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Abstract 

We analyze in detail the size of next-to-leading order corrections to hard exclusive meson 
production within the collinear factorization approach. Corrections to the cross section 
are found to be huge at small xb and substantial in typical fixed-target kinematics. With 
the models we take for nucleon helicity-flip distributions, the transverse target polarization 
asymmetry in vector meson production is strongly affected by radiative corrections, except 
at large xb- Its overall size is very small for p production but can be large in the u channel. 
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1 Introduction 



Generalized parton distributions (GPDs) have developed into a versatile tool to quantify important 
aspects of hadron structure in QCD. In particular they contain unique information on the transverse 
spatial distribution of partons [1] and on spin-orbit effects and orbital angular momentum inside 
the nucleon (2j [3] . Deeply virtual Compton scattering is widely recognized as the process providing 
the theoretically cleanest access to GPDs, with a wealth of observables calculable in the large Q 2 
limit [4] and with the calculation of the hard-scattering subprocess now pushed to next-to-next-to- 
leading order (NNLO) accuracy in a s [5\. A quantitative theoretical description of exclusive meson 
production remains a challenge. It would offer the possibility to obtain important complementary 
information, difficult to obtain from Compton scattering alone. Perhaps most importantly, vector 
meson production is directly sensitive to gluon distributions, which in the Compton process are a s 
suppressed relative to quark distributions and only accessible through scaling violation (just as in 
the well-known case of inclusive deep inelastic scattering). Given in addition the large number of 
channels that can be studied and the wealth of high-quality data in a wide range of kinematics from 
collider to fixed-target energies 0(7], it should be worthwhile to try and push the theory description 
of exclusive meson production as far as possible. 

In this work we study the exclusive production of light mesons at large photon virtuality Q 2 
within the framework of collinear factorization [8j. In Bjorken kinematics, the process amplitude can 
be approximated by the convolution of hard-scattering kernels with generalized parton distributions 
and the quark-antiquark distribution amplitude of the produced meson. The hard-scattering kernels 
have been calculated to 0(a 2 ), i.e. to next-to-leading order (NLO) accuracy [9], [10l [TT] . The aim of 
the present paper is to investigate in some detail the size of the NLO corrections compared with the 
leading-order (LO) results, on which phenomenological studies have so far relied. 

The collinear factorization approach provides an approximation of the leading helicity amplitudes 
for meson production in the Bjorken limit, up to relative corrections of order 1/Q 2 . These power 
corrections cannot be calculated systematically (and in fact the derivation [8] of the factorization 
theorem suggests that these corrections do not all factorize into hard-scattering kernels and nonper- 
turbative quantities pertaining to either the nucleon or the produced meson). One particular source 
of power corrections can however readily be identified, namely the effect of the transverse momentum 
of partons entering the hard-scattering subprocess, which in the collinear approximation is neglected 
in the calculation of the hard-scattering kernel. A number of approaches include these kx effects, 
in particular the studies in |12} [T3] based on the modified hard-scattering formalism of Sterman et 
al. |14j . and calculations like [15] which are based on the color dipole formulation. In the work by 
Martin, Ryskin and Teubner [16J, parton-hadron duality is used to model the meson formation and 
thus the transverse momentum of the hadronizing quarks is included in the calculation, whereas the 
transverse momentum of gluons in the proton is treated within high-energy kx factorization. The 
studies just quoted agree in that transverse momentum effects result in substantial power corrections 
to the collinear approximation for Q 2 up to several GeV 2 . Unfortunately, the calculation of full 
NLO corrections in a s remains not only a practical but also a conceptual challenge in all of these 
approaches, so that the perturbative stability of their results cannot be investigated at present. (The 
approach of Sterman et al. takes partial account of radiative corrections, resumming a certain class 
of them into Sudakov form factors.) 

A consistent simultaneous treatment of radiative and power corrections being out of reach at 
this time, a possible strategy is to study the NLO corrections in the collinear approximation and in 
particular to identify kinematical regions where these corrections are moderate or small. There one 
can then use with greater confidence formulations incorporating power corrections. In this spirit the 
present investigation should be understood. We will study both the cross section for meson production 



3 



from an unpolarized target and the transverse target polarization asymmetry. This asymmetry is one 
of the few observables sensitive to the nucleon helicity-flip distributions (in particular for gluons) and 
hence to the spin-orbit and orbital angular momentum effects mentioned above. We will in particular 
see whether corrections tend to cancel in this polarization asymmetry, as is often assumed. 

In the bulk of this paper we concentrate on the production of vector mesons. In Sect. [2] we set up 
our notation and recall important properties of the hard-scattering kernels at NLO, as well as giving 
a one- variable representation of these kernels after Gegenbauer expansion of the meson distribution 
amplitude. In Sect. [3] we specify the model of the generalized parton distributions H and E we use for 
our numerical studies. The size of radiative corrections involving convolutions with distributions H 
is then studied in Sects. [Hand [5] for small and large xb, respectively, and the convolutions involving 
distributions E are quantified in Sect.[6j In Sect.[7]we then look at the NLO corrections at the level of 
the observable cross section and polarization asymmetry. A brief study of exclusive pion production 
in Sect. [8] complements our work, and in Sect. [9] we summarize our main findings. A number of more 
lengthy formulae is collected in appendices. 



2 Hard-scattering kernels 

In the main part of this paper we are concerned with exclusive production of a vector meson 

l*(q)+p(p)^V(q')+ P (p') (1) 

in the limit of large Q 2 = —q 2 at fixed Bjorken variable xb = Q 2 /{2p • q) and fixed t = (p — p') 2 . To 
leading order in 1/Q, the amplitude for longitudinal polarization of photon and meson can be written 
as 
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with z = 1 — z, N c = 3, and the electromagnetic fine structure constant a. Throughout this paper 
we work with nj = 3 active quark flavors. The proton matrix elements F are parameterized by 
generalized parton distributions, 
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for quarks and gluons, where we use the conventions of [T7J. Here n is a light-like auxiliary vector, 
£ = #b/(2 — xb) is the skewness variable, and m p denotes the nucleon mass. We have further 
introduced the combination 



F*+Hx,t,t)=F*(x,Z,t)-F*(-x,Z,t) 



(4) 



with positive charge conjugation parity. In ^) we have arranged the terms containing quark distri- 
butions into the flavor singlet 



F S = F u(+) + pd{+) + F s(+) 



(5) 



4 



V 




n 



Figure 1: Example graphs for the hard-scattering kernels T a , and T„ at order a 2 . 
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of the p, uj and <j). The meson distribution amplitudes 4>v( z ) are normalized as J dz<j)y(z) = 1, and 
the decay constants have the values f p = 209 MeV, f u = 187 MeV, = 221 MeV [18]. We finally 
have hard-scattering kernels in ([2]) , where T g goes with gluon and T a , go with quark distributions 
in the proton. In the graphs for T a quark lines connect the proton and meson side, whereas in the 
graphs for Tj, the proton and meson side are only connected by gluon lines. T& thus starts at order a 2 s 
and only goes with the quark singlet distribution F s . Example graphs for the three kernels at NLO 
are shown in Fig. [U We will refer to T g , T a , as the gluon, the quark non-singlet, and the pure 
quark singlet kernel, respectively. 

For better legibility we have not displayed the dependence on the renormalization and factorization 
scales in ([2]) . The renormalization scale [i R appears as argument of a s and through explicit logarithms 
in the hard-scattering kernels T. The kernels further contain logarithms of the respective factorization 
scales [Ida and hgpd for the meson distribution amplitude and the generalized parton distributions. 
The NLO kernels in [1Q(, lllj are given for a common factorization scale fip = hdA = Hgpd- We can 
restore the individual logarithms of hda and hgpd from the requirement that within the calculated 
precision the process amplitude ([2]) must be independent of these scales. As an example consider the 
term 

d f 1 

— r— / dz<pv(z;^DA)T a (z,X,^]a a (fi R ),fi R ,HGPD,tiDA,Q) 
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0) 



where the scale dependence of <\>y{z\ Hda) is given by the ERBL evolution equation [19J. At lead- 
ing order this gives a term d / d(ki fJ, DA ) (f>v(z; [Ada) of order a s , whose convolution with the 0(a s ) 
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part of T a must cancel against the contribution from explicit logarithms of hda in the 0(a 2 ) part 
of T a . An analogous argument holds for the dependence on ^gpd, with the complication that 
the gluon and quark singlet distributions mix under evolution. More precisely, the convolution of 
d/d(ln/i^p D ) F s (x,^,t; hgpd) with the 0(a s ) part of T a cancels at 0{a 2 s ) against the contributions 
from logarithms of [j-gpd in T a and in T g . Likewise, the convolution of d/d(ln n GPD ) F 9 (x, £, t; hgpd) 
with the Born term of T g cancels at 0(a 2 ) against the contributions from logarithms of hgpd in T g 
and in the pure singlet kernel T b . We have explicitly checked that the scale dependence of the hard- 
scattering kernels given in [11] cancels in the process amplitude ([2]) as just described, using the LO 
evolution equations for GPDs given in App. O 

Separating the [id A and hgpd dependence, we can write the kernels as 
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where y = 1 — y and we use the standard notation 
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The functions IC g , K. b and /C a are independent of Q 2 and the renormalization and factorization scales. 
They contain factors Cp or C A but not (3q. Their expressions can be found in [11 j . taking into account 
that the kernels T g and T b here are denoted by T g and T( + ) there, and that 
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Note that the pure singlet kernel T b does not contain logarithms of hda and hr at 0(a 2 ), since there 
is no Born level contribution against which they could cancel in the scale dependence of the process 
amplitude. There is however a logarithm of hgpd, since the corresponding derivative of the Born level 
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convolution of T g with F 9 contains a term going with the quark singlet distribution F , as already 
mentioned after ([9]). 

The kernels in (I10p have singularities for real- valued arguments. One readily finds that x/£ = 
(s — u)/Q 2 , where s and u are the Mandelstam variables for the parton- level subprocess 7*5 — > (qq)q 
or 7*g — > (qq)g. The prescriptions s + ie for the s-channel and u + ie for the u-channel singularities 
thus instruct us to take x + ie for x > and x — ie for x < 0. Correspondingly, the second argument 
(£ — x)/(2£) of T 5 , Tfe and T a must be taken with — ie for x > and +ie for x < 0. In T a (z, — x, £) the 
second argument of T a is (£ + x)/(2£), which has to be taken with — ie for x < and +ie for x > 0. 
We remark that, as it is written, the ie prescription in [IT] for the gluon and the pure singlet kernel 
is correct for x > but incorrect for x < 0. Likewise, the prescription given in [10} [11] for the quark 
non-singlet kernel is correct for x > if the corresponding argument is (£ — x)/(2£) and for x < if 
the argument is (£ + x)/(2£), but incorrect in the other caseslj 



2.1 Gegenbauer expansion 

Let us expand the meson distribution amplitude on Gegenbauer polynomials, 

00 

M*;/0 =6z(l-z)Y,a n ^)C^ /2 (2z-l), (15) 
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where do = 1 according to the normalization condition J dz <\>y{z) = 1. To leading order, the 
Gegenbauer coefficients evolve as 
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where a s (fi) is the running coupling at one-loop accuracy. One has 7„ ~ 4Cp ln(n+ 1) within at most 
6% for all n. For V = p,u,(f) only coefficients a n with even n are nonzero due to charge conjugation 
invariance, and in all subsequent expressions of this paper we consider n to be even. Calculations of 
the distribution amplitudes in models or on the lattice typically give values for the first or the first 
two nonvanishing moments, see e.g. |2U[ I2"l| I22| . so that a truncated version of the expansion (|15p is 
very often used in phenomenological studies. Convolution with individual terms in (|15p also allows 
us to reduce the hard-scattering kernels for meson production to functions of a single longitudinal 
variable. More precisely, we can rewrite the process amplitude ([2]) as 
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with convolutions in x 
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1 We thank Dima Ivanov for discussions on this point. The numerical results in were obtained with the correct 
prescription. 
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which depend on £ and i, and logarithmically on Q 2 and on the factorization and renormalization 
scales. At order c? s the dependence on /ir and on hda cancels in each separate convolution, while the 
dependence on hgpd cancels in J 7 ^ and JF® and in the sum + + as discussed after 

Q. In analogy to (|19p we define convolutions TL and £ for the individual distributions H and E in 
(J3]). The kernels T Si „, T Qjn , Tfe jn are obtained from T g , T a , Tb by multiplying with 6z(l — z) C^ 2 (2z — 1) 
and integrating over z. For n = we find 



2£ 



(£, — x — ie)(l; + x — ie) 



2£ 



T 0j n(x,£) = -3CFO.S 



£ — X 



£ — x — ie 



1 ~\~ . t a n 

4tt 



S, — X 



(20) 



with 

*s,o(y) 



2C A (y 2 + y 2 )-C F y 



ln2/ ln^ + ^ln 



1 



+ ( - + l-4y ] lny 



V In y 
2 y 



7T 



2 (y - V) Li 2 y - 4yy ( 3 Li 3 y - In y Li 2 y - — In y 



*6,o(y) = 2(y - y) ^ 
y 



8y In y + 



y 



In 



Q' 2 



^GPD 



2y ) In 2 y + 2(y - y) Li 2 y 
ln 2 y 



+ {y -> y} , 



+ (y - y) — ^- + 4 Li 2 y - {y -> y} 

y 



ia,o(y) = A 



19 . . Q 21 

6 



(3 + 21ny) In 



<3 2 



77 



^GPD ® 



1 



3 In y + In 2 y 



(2C F - C A ) | - i - 4(2 - 3y) In y + 2(1 - 6y) In y + 4(1 - 3y) (Li 2 y - Li 2 y) 

r _ vr 2 _ 

+ 2(1 - 6yy) 3(Li 3 y + Li 3 y) - lny Li 2 y - lny Li 2 y - — (lny + lny) 



(21) 



The corresponding kernels for n = 2 and n = 4 are given in App. [Bj The ie prescription to be used 
in f[20|) is the same as specified at the end of the previous subsection. This implies that in t gtTl (y), 
tb,n{y)i t a ,n(y) and t a ,n(y) one has to take ln(y — ie), Li 2 (y + ie) and Li 3 (y + ie) for y < 0. For the 
gluon and pure singlet kernel, which dominate in process amplitudes at small £, we have in particular 
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in the region y < 0. In the limit y —* all three expressions in (|2ip contain singular terms proportional 
" " L ,! : /f 19j) we should however consider {yy)~ 1 t g>n (y), y^ 1 t a>n (y) and 
20]) . With the appropriate prescription, these kernels contain terms which 
_1 ln m (y — ie), where m = 0, 1, 2. 



to lny and ln^ ?/. For the convolution 
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3 Model for the unpolarized GPDs 

It is difficult to study the impact of NLO corrections at the level of the hard-scattering kernels given 
in the previous subsection, especially since they are not smooth functions but distributions with 
singularities at y = 0. We will therefore use model GPDs to investigate the radiative corrections at 
the level of the convolution integrals (|19j) . The aim of this work is not a systematic improvement of 
existing models, nor a detailed exploration of model uncertainties on observables in exclusive meson 
production. We do however require that the models we use are consistent with known theoretical 
requirements and basic phenomenological constraints. 

For H q and H 9 we adopt the widely used ansatz of [261 EZj based on double distributions, where 
a £ dependence is generated according to 

P-$a) /i (2) (/3,a) H q ( + \(3,0,t) , 
-ta)h {2) ([3,a)H 9 ((3,0,t) (23) 

T(26 + 2) [(l-|/?|) 2 -a 2 ] b 

2 2&+i r 2(6 + l) (1-^1)26+1 • ^> 

taken as 

H^{x,0,t) = q v {x)exp[tf qv (x)] +2q{x)exp[tf q (x)) , 

H°(x,0,t) =xg(x)exp[tf g (x)] (25) 

for x > 0, with the values for x < following from the symmetry properties of the distributions. Here 
q v (x) = q(x) — q(x), q(x) and g(x) are the usual unpolarized densities for valence quarks, antiquarks 
and gluons, for which we take the CTEQ6M parameterization [29]. This parameterization has an 
identical strange and antistrange sea, so that s v {x) = 0. The ansatz (j23[) is taken at a starting scale 
jUo and then evolved with the LO evolution equations given in App. [Cj For the studies in Sects. 2] 
and [5] we take fio = 1.3 GeV, which is the starting scale of evolution for the CTEQ6M densities. In 
Sects. [6] and [7] we will instead take fj,Q = 2 GeV, since this will allow us to use the results for the t 
dependence of valence distributions obtained in [28J. 

For the t dependence in the ansatz (I23p we follow the modeling strategy of [27j and take an 
exponential behavior in t with an x dependent slope. For valence quarks we take the slope functions 

f qv (x) = a' v (l- xf In - + B Qv (1 - x) 3 + A qv x(l - x) 2 (26) 
x 

with parameters a' v = 0.9 GeV -2 and 



H^\x,^,t)= I dp I 
J-i J-i- 



H 9 (x,tt) 



dp 
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da 8(x 
da 8(x 



with 

h {b) (P,a) = 

The distributions at zero skewness are 
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A Uv = 1.26 GeV" 2 , B Uv = 0.59 GeV~ 2 , 

A dv = 3.82 GeV" 2 , B dv = 0.32 GeV" 2 , (27) 



from [28]. We recall the sum rule 

= J dxH q {x,0,t) = J dxq v (x)exp[tf qv (x)] , (28) 

from which one obtains the electromagnetic Dirac form factors of proton and neutron by appropriate 
quark flavor combinations. Together with the CTEQ6M distributions at (mq = 2 GeV, the ansatz in 
(|26p and (|27p gives a good description of the data for these form factors. For gluons we take a slightly 
simpler form than (126j) and set 

f g (x) = a'Jl - xf In - + B g (l - xf . (29) 

For the parameters we take 

a' g = 0.164 GeV -2 , B g = 1.2 GeV" 2 (30) 

so as to match recent HI data on J/* photoproduction, whose t dependence is well fitted by [30] 



da 

— oc exp 



6 + 4^1n^ 



(31) 



with central values b = 4.63 GeV -2 and a' g = 0.164 GeV" 2 for W = 90 GeV. To connect §T$) 
with (f29|) we have used the approximate relation da/dt oc \H 9 (£, t)| 2 , which is obtained when only 
keeping the imaginary part of the tree-level amplitude, where 2£ = (Mjj^,/W lv ) 2 in terms of the 
7£> cm. energy. With the ansatz (|23p one approximately has H 9 (^,^,t) oc exp 

[tf g (2£)] for the t 

dependence of the GPD [13]. 

Whereas information on valence quark GPDs can be obtained from the sum rules (j28|) and infor- 
mation on gluon GPDs from J/^f production, almost nothing is so far known about the t dependence 
of GPDs for antiquarks. As a simple ansatz we shall take their slope functions equal to those in the 
valence sector, 

fu = fu v 5 fd = fd v j fs = fd v j (32) 

bearing in mind that it remains an outstanding task to develop more realistic models. 



3.1 Nucleon helicity-flip distributions 

The nucleon helicity-flip distributions E q and E 9 are less-well known than their counterparts H q and 
H 9 , because their values at £ = and t = cannot be measured in inclusive processes and are thus 
subject to considerable uncertainty. 

The model described in this subsection refers to a scale of /1q = 2 GeV. We make a double 
distribution based ansatz 

P-£a) hW(J3,a) (/?, 0, t) , 
(3-£a)h i2 \[3,a)E 9 {(3,0,t) (33) 



E^+\x,tt) 
E 9 (x,£,t) 



dp 



dp 



1+1/31 



da 8(x 



1+1/31 



da 8(x 
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as in (f23l) . and for x > set 

(x, 0, t) = e g „ (a?) exp [t 5 ,„ (x)] + 2e 5 (s) exp [t^(x)] 



E g (x, 0, i) = xe g (x) exp [t<7 g 



x 



(34) 



with the corresponding values for x < determined by the symmetry properties of the distributions. 
For the forward limit of the valence distribution we take 



whose normalization factor 



ensures the sum rules 



e qv {x) =K q N(a v ,(3 qv )x~^(l-x)^ , 
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(35) 
(36) 

(37) 



where k u ~ 1.67 and K d ~ —2.03 are the contributions of u and d quarks to the anomalous magnetic 
moment of the proton. For the functions controlling the t dependence we take the same form as 
in d26D, 

g qv (x) = o/(l - x) 3 In - + D qv (1 - x) 3 + C Qv x(l - xf . 



(38) 



With the parameters a v = 0.55, a' v = 0.9 GeV 2 and 



f^U — 3.99 . 

P d = 5.59 ; 



C Uv = 1.22 GeV" 
C dv = 2.59 GeV" 



D Uv = 0.38 GeV , 
D dv = -0.75 GeV" 2 



(39) 



from [28j one obtains a good fit to the electromagnetic Pauli form factors of proton and neutron via 
the generalization of the sum rule (|37[) to finite t. 

For the forward limit of the distributions of antiquarks and gluons we make the same simple ansatz 
as in (f35|) . 



eg(x) = k q x- a i(l - x)h , 
and for the t dependence in the gluon sector we set 

9g{x 



e g {x) = kgx~ as {l - xf<> , 



a' g {l-x) 2 ln-- + D g (l-x) 2 , 



(40) 



(41) 



in analogy to the form (|29p we used for H 9 . We presently have not no phenomenological information 
on these distributions, but two theoretical constraints. There is a condition that ensures positive 
semidefinite densities of partons in the transverse plane [31] . which with our ansatz for the GPDs 
reads [28] 



"eg(x)~ 


2 

< 8errip 




q(x) _ 




~e g (x)~ 


2 

< 8em 2 




. S(x) _ 





[fq(x) - gg(x)] , 
[fg(x) -9g{x)] 



(42) 
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if we neglect for simplicity the polarized antiquark and gluon distributions compared with the unpo- 
larized ones. On the other hand we have the sum rule 



0=1 dxE 9 (x, 0,0) + V f dxxE q (x, 0,0) 
Jo q J-i 

= [ dxxe 9 (x) + y2f dx x[e qv (x) + 2e q (x)] (43) 

JO q 

following from the conservation of the energy-momentum tensor. For the parameters in (|4ip we take 

a' g = 0.164 GeV" 2 , D g = 1.08 GeV~ 2 , (44) 

with a' g as in (|30p and D g slightly smaller than its counterpart B g for H 9 , so that the positivity 
condition f|42[) can be fulfilled. Assuming a similar small- x behavior of the distributions for proton 
helicity-flip and non-flip, we take in (140p the values a q = 1.25 and a g = 1.10, which we obtain when 
fitting the CTEQ6M distributions to a power law in the x range from 10" 4 to 10" 3 . 

Since it turns out that the transverse target polarization asymmetry in p production is very 
sensitive to the details of the helicity-flip distributions, we will explore two model scenarios in our 
numerical studies: 

1. a scenario where the sea quark distributions e q behave similarly to the valence distributions e qv . 
For the t dependence we then take gu(x) = g Uv {x) and gs{x) = gd v {x). The parameters k q in 
(|40p are taken such that second moments at t = fulfill 

fpdxxeg(x) = ^dxxg{x) 
J dx xe qv (x) dx xq v (x) 

for q = u,d, where the ratio on the r.h.s. is taken from the CTEQ6M parameterization at fi = 
2 GeV. Its value is 0.095 for u and 0.30 for d quarks. This fixes the values of k q N _1 (a q — 1, f3g) 
with N given in ([36"p . For the strange distribution we set e s = e s = 0, and k g N~ 1 (a g — l,/3 5 ) 
is then fixed by the sum rule (|43p . 

The powers (5g and /3 g controlling the large-x behavior are finally taken to have the smallest 
values for which the positivity condition (|42p holds in the range x < 0.9 (for higher x even the 
unpolarized densities are so uncertain that we do not insist on the positivity conditions to be 
fulfilled). 

2. a scenario where e q behaves similarly to the gluon distribution e g . The t dependence is now 
modeled by taking g q {x) = g g {x) for q = u,d, s. For the second moments we impose 

Jldxxe q {x) _ J^dxxqjx) 
Jq 1 dx xe g (x) Jq 1 dx xg(x) 

for the three light quark flavors, where with the CTEQ6M distributions the r.h.s. is equal to 
0.064, 0.083, 0.036 for u, d, s, respectively. We now have a nonzero e s = e s . The values of 
k q N~ 1 (a q — 1, (5 q ) and k g N~ 1 (a g — 1, /3 g ) are taken to fulfill both (|4"6"]1 and ([4^1 . and the powers 
j3q, (5 g are set to the minimal values for which positivity holds in the range x < 0.9. 

The parameters resulting from this modeling procedure are collected in Table [H and the distributions 
at £ = and t = for model 1 are shown in Fig. [2j 
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Table 1: Parameters in the ansatz (I40p for different parton species a in the two models described in 
the text. The values for valence quarks apply to both models, with normalization parameters given 
by kg v = K q N(a v , (3g v ) according to (|35|) . The last line gives the second Mellin moment at \i = 2GeV 
in the forward limit. 
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Figure 2: The forward limits e a (x) of the nucleon helicity-flip distributions at /i = 2GeV for different 
parton species a in model 1. 

We find that in model 2, both sea quark and gluon distributions are nearly zero (so that we do 
not attach importance to the unrealistically small value of (3 g obtained with our above procedure). 
Their smallness can be traced back to the small value of the flavor singlet integral 

f dx x [e Uv (x) + e dv (x)] = 0.008 (47) 
Jo 

in the valence sector of our ansatz. In model 2, the distributions eq and e g have the same sign as a 
consequence of (146p and due to the sum rule (|43p can only be tiny. Somewhat larger distributions for 
sea quarks and gluons are obtained in model 1, where they have opposite sign because of (|45p . 

We note that the parameters (|39p we have taken for the valence part of E q are by no means 
precisely determined by a fit to the Pauli form factors: alternative fits in [28] gave a similarly good 
description of the form factor data, with some variation of the resulting value of the integral in (|47p . 
Nevertheless, any model where e Uv and ed v have similar shapes and no zeroes in x will yield rather 
small values of this integral, given the strong cancellation between u and d quark contributions in the 
moment dx [e Uv (x) + ed v (x)] = k u + « —0.36. It would be interesting to explore how much the 
integral ([17J) and as a consequence the sea quark and gluon distributions can vary in realistic models, 
but such an investigation is beyond the scope of this work. 
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We end this section by quoting the values for the total angular momentum carried by quarks and 
antiquarks of a given flavor in our model, given by 

1 f 1 

J q = - dxx[H q (x,Q,0) + £ 9 (x,0,0)] (48) 

2 7-1 

according to Ji's sum rule [2j. With the parameters in Tableland the CTEQ6M distributions we 
find 

J u = 0.25 , J d = -0.01 , (model 1) 

J u = 0.24 , J d = 0.03 (model 2) (49) 

at the scale \x = 2 GeV of our model. We note that this is in rather good agreement with the 
results of recent lattice calculations, with J u = 0.214(16) and Jd = —0.001(16) reported in [32], and 
J u = 0.33(2) and J d = —0.02(2) in [33]. Let us reiterate that with just two sets of model parameters 
we cannot exhaust the range of possible scenarios but only provide two representatives that are 
consistent with presently known constraints. As just discussed, the relative smallness of sea quark 
and gluon distributions compared with the nucleon helicity conserving case should however be typical 
of a rather wide class of models. 

4 Vector meson production at small xb 

We now study numerically the importance of NLO corrections in vector meson production. Here 
and in the following sections we use the two-loop strong coupling for nt = 3 flavors with a QCD 
scale parameter A® = 226 MeV. This value corresponds to = 326 MeV, A^ 5 ) = 372 MeV and 
to a^(Mz) = 0.118 when matching at m c = 1.3 GeV and mj = 4.5 GeV, which are the values 
used in the CTEQ6M parton analysis [29|. We also take nj = 3 fixed in the evolution and the 
hard-scattering kernels. Taking rif = 4 with massless charm or rif = 5 with massless charm and 
bottom would not be a good approximation for the rather moderate values of Q 2 we will discuss 
for fixed-target kinematics. On the other hand, taking rif = 3 and neglecting charm altogether is 
admittedly not a good approximation for the larger Q 2 relevant in collider kinematics. However, with 
c4 3 ) = 0.164 compared to = 0.178 at /J, = 10 GeV we expect that this inaccuracy will not affect 
the conclusions at high Q 2 we shall draw from our studies. 

We have performed the evolution of the GPDs at LO using the momentum-space evolution code 
of [33]. As explained in Sect. [2] taking LO evolution together with the NLO hard-scattering kernels 
is sufficient to obtain scale independence of the process amplitude up to uncalculated corrections of 
order a^. With the input scale of evolution not taken too small, NLO evolution effects should be 
rather moderate at the Q 2 values relevant in fixed-target kinematics, whereas our general conclusions 
for high Q 2 and small xb will again not depend on this level of detail. We note that the NLO kernels in 
momentum space are available in the literature [35] , but their considerable length makes it difficult to 
implement them in a fast numerical evaluation. For including NLO effects in the evolution it should 
be more efficient to use the Mellin space approach recently followed for deeply virtual Compton 
scattering in [5]. 

Here and in the following section we consider the convolutions of hard-scattering kernels with 
GPDs at t = 0. For nonzero £ = x#/(2 — xb) this should be understood in the sense of an analytic 
continuation, since the physical region for meson production is —t > Am 2 ^ 2 /(l — £ 2 ) in Bjorken 
kinematics. To explore the importance of NLO corrections we do not see this as a shortcoming. 

Let us start our discussion with the gluon and quark singlet sector. Here and in following we shall 
always present the convolutions (119D for Gegenbauer index n = unless indicated otherwise. In Fig. [3] 
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Figure 3: LO and NLO terms of the convolutions in the gluon and quark singlet sector at Q = 4 GeV. 
The scales are set to fj,R = hgpd = [J-DA = Q- The NLO terms are for Gegenbauer index n = unless 
specified explicitly. Here and in the following plots the label "NLO" denotes the 0(a 2 s ) part of the 
convolutions, whereas the sum of 0(a s ) and 0(a 2 s ) terms is labeled by "LO+NLO". 
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Figure 4: LO terms and the sum of LO and NLO terms of the convolutions in the quark non-sing 
sector at Q = 4GeV, with = ^gpd = Hda = Q- 
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Figure 5: Dependence on the common scale fi = [ir = hgpd for the sum of convolutions in the gluon 
and quark singlet sector. 
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Figure 6: Cross section for j*p — > pp with longitudinal photon polarization. Bands correspond to 
the range Q/2 < \i < 2Q and solid lines to p = Q. We also show the power-law behavior a oc l^ - 88 
(with arbitrary normalization) obtained from a fit to data in the range 0.001 <x^^ 0.005 |36j . 

we show the LO and NLO pieces of the convolutions for the scale choice /xr = pgpd = /ADA = Q- The 
size of corrections at small xb is dramatic: we have large NLO corrections with opposite sign compared 
to the LO term for 7i 9 , and a similarly large NLO contribution from Tl s ^ with sign opposite to the 
LO result for TL s ^ a \ In the sum of gluon and quark singlet terms, the NLO corrections drastically 
reduce the LO result or even lead to a change of sign between LO and the sum of LO and NLO 
results. We also observe that for higher Gegenbauer index the NLO corrections tend to be even more 
important. Note that the LO term of the convolutions is the same for all n as can be seen from 
(|19p and (|20p . The size of NLO corrections in 7i s ^ is comparatively moderate, at least for lower 
Gegenbauer moments. The same is seen for the quark non-singlet convolutions in Fig. SJ Of course, 
the gluon and quark singlet terms will dominate meson production at small in those channels 
where it is allowed by the meson quantum numbers. 

In Fig. [5] we explore the influence of the scale choice by varying p^ = pgpd simultaneously. For 
xb = 2 x 10 -3 we find an indication for the onset of perturbative stability at Q = 7GeV but not yet 
at Q = 4GeV. For xb = 2 x 10~ 2 the situation is less severe, with moderate corrections in a wide 
p range already at Q = 4GeV. In contrast, when going down to xb = 2 x 10 -4 we find very large 
corrections even at Q = 7GeV. We have checked that the conclusions in the respective kinematics 
do not change when we vary pgpd while keeping Pr = Q fixed. 

Figure [6] shows how the perturbative instability we observed in the convolutions affects the longitu- 
dinal cross section for p production. Here we have taken the asymptotic form of the meson distribution 
amplitude, i.e. set a n = for n > 2. In the NLO result for the cross section we have squared the co- 
herent sum of LO and NLO terms in the process amplitudejl i.e. we have taken |A4lo + A4nlo| 2 - We 
see that the NLO corrections severely decrease the LO result. As a consequence of the cancellations 
between LO and NLO contributions, the scale dependence of the cross section does not decrease. We 
also show in the figure the power-law behavior a oc iy - 88 obtained from a fit to data in the range 
0.001 <x B < 0.005 [36]. As observed in [13], a double distribution model with the CTEQ6M distri- 
butions as input lead to a rather good description of this experimentally observed energy dependence 

2 We thus keep terms of O(a^) in the cross section, although the accuracy of the NLO calculation is only up to 0(o? s ). 
This should not be seen as a problem, as it will not make a considerable difference in situations where perturbative 
corrections are moderate, whereas in situations where NLO corrections are huge we would neither trust the cross section 
with or without the partially included 0(o? s ) terms. 
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Figure 7: The factorization scale dependent and independent terms of Im t g ^ n and Im t bjU as specified 
in (|54p . shown for n = and y < 0. 



if the cross section is evaluated at LO. With the strong cancellations from the 0(a 2 s ) corrections, one 
obtains an NLO result whose energy behavior is much too weak. 

Let us discuss how the huge size of corrections can be understood at an analytical level, following 
the line of argument given in [37] . Using (|8ip and (|82p we can approximate the hard-scattering 
kernels for large negative y as 



- lml g (z,y) = AC a 



7T 



2 1 



m(zz) + In — ^ — 
^GPD 



-lml b (z,y) = 4 

7T 



1 — m(zz) 



In 



Q 2 



t^GPD 



y + o(i), 

+ 0(y- 1 ) 



Rel g (z,y)=0(l), 
Rel b (z,y)=0(y' 1 ), (50) 



where here and in the following the order of corrections is given up to powers of lny. The quark 
non-singlet kernel is subleading compared with the pure singlet one, 



0(y 



(51) 



where we have divided X a (z, y) by y corresponding to the prefactor in the complete kernel (|10p . From 
(|50|) we readily obtain 



1 



Im t g , n (y) = -4C A 



TT 



-lmt bn (y) = 4 

7T 



hi 



Q 2 



^GPD 



Cn + 1 — In 



Q 2 



^GPD- 



y + o(i) 

+ 0(y~ l ) 



with constants 
c = 2, 



11 

C2 = y 



„ 137 _ 

3.7 , C4 = « 4.6 , 

30 



/ cfe ln(zz)C^ 2 (2 
Jo 



z-1) 



(52) 



(53) 



that increase with the Gegenbauer index n. In Fig. [7] we show for the case n = that these approxi- 
mations become very good for increasing \y\, where we have decomposed the exact kernels as 



Q 2 



Im * 9 ,n(y) = lmt g.n(y) + lmt g.n(y) ln ~2 

^GPD 
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Im V(y) = lmt b, n (y"> + lmt b, n (y) lr W 

^GPD 

Let us now rewrite the convolutions of kernels and GPDs in terms of the variable lo = x/£, 



G 



Q 1 



(54) 



Im 



H 9 n + H s ® 
2 



-6a., 



a* 



nHS(£,£,t) + - ± Im / dw 



1 — to — ie 



l-u>-ie \ gW) „ , 



1+co 



1-LO 



IE 



(55) 



For lo > ujq with some wo > 1 we can use the approximation (|52p of the hard-scattering kernels, and 
further approximate 1 + lo ~ to in the first term on the second line. This gives 



1 



67ra s 



7T 



Im 



n 9 +n^ b) 



TT 



dco ... 



da; < Ca 



In 



Q' 2 



^GPD 



H 9 (^,£,t) 



LO 



+ c F 



c n + 1 - In 



Q 2 



^GPD 



(56) 



where the integral over u; on the first line is to be taken with the unapproximated integrand from 
(|55|) . It grows with £ like H 9 (to^,^,t) or £H S (Lu£,£,t) but lacks the enhancement due to the upper 
limit l/£ of the integral on the second line. Restricting our discussion to t = for simplicity, we can 
for sufficiently large lo neglect the effect of skewness in the GPDs and then have 



H s (^t^)^S{Lod)=Y J [9K)+?K)] 



(57) 



where S{x) is the usual quark singlet distribution. In a very rough approximation one may treat 
xg{x) and xS(x) as constant at small x. In (|56p one then has loop integrals J cIlo/lo for both the 
gluon and the quark term, which generate large logarithms ln(u;o£) for l/£ S> wq- These logarithms 
are of BFKL type and correspond to graphs with i-channel gluon exchange in the hard-scattering 
kernel, such as those for and T g in Fig. [TJ 

In a phenomenologically more realistic approximation one has xg(x) ~ ax~ x at small x and a 
similar behavior with different values of a and A for xS(x). This gives 



or 



do; a;' 



-A-l 



(58) 



for l/£ 3> wo, when the bulk of the integral comes from the region where the small-x approximation 
of the gluon density is valid. With A being rather small for the gluon distribution in a wide range of 
the factorization scale, the term ([58]) has the same power behavior £~ A as the Born term H 9 (£,£, 0) 
in (|56p but is numerically enhanced by 1/A. A contribution analogous to (|58p is obtained from the 
quark singlet term in (|56p and comes with a similar enhancement. 

Concerning the choice of factorization scale, it is clear that the size of the corrections in (|56p is 
decreased if hgpd is taken smaller than Q. It is also clear that no scale choice can eliminate both 
the gluon and quark singlet contribution in this expression. To make at least the gluon term for 
n = disappear one needs Hqpd = e~ 2 Q 2 m 0.14 Q 2 . For a wide range of Q 2 this is outside the 
perturbative region or at least so low that the quark singlet distribution has a rather small power A 
and can thus give important corrections. We note that previous analyses of vector meson production 
at small xb have argued for a factorization scale well below Q 2 , based on different estimates of the 
typical virtualities in the leading-order graphs |151 138] . We also note that the fiR dependent term 
Po \n{^i 2 R / ijlq PD ) in the gluon kernel (fTT|) does not appear in the approximation (f50l) which dominates 
the convolutions at small sg. The choice of can thus not cure the huge NLO corrections we have 
discussed. 
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5 Vector meson production at moderate to large xb 



Let us now investigate the NLO corrections in typical fixed-target kinematics, as it is accessible 
at HERMES, JLab and COMPASS. We take again t = and for definiteness present estimates at 
Q 2 = 4GeV 2 . For larger Q 2 , which will in particular be accessible with the JLab energy upgrade to 
12 GeV, the corrections are in general smaller. 

In Fig. [8] we compare the LO and NLO parts of the convolution integrals. In the gluon sector 
we find no simple picture, with relative corrections that are typically moderate but become large 
for JieTC 9 at smaller xb and for ImTC 9 at larger xb- For the quark singlet the situation is similar 
to the one in the small-xs region, i.e. we have rather large NLO corrections from 7i s ^ with sign 
opposite to the LO part of 7i s ^ , whereas the NLO corrections in 7i s ^ are smaller. Adding gluon 
and quark singlet contributions, we find that for n = the NLO corrections are of reasonable size for 
the imaginary part. For the real part at lower xb, the corrections are however large and of opposite 
sign compared to the Born term. We note that the convolutions Ti satisfy a dispersion relation in 
1/xb for fixed Q 2 and t [39]. In this representation their real parts at a given xb are sensitive to the 
imaginary part at smaller values of xb, where the NLO corrections rapidly increase as we have seen 
in the previous section. Turning to the quark non-singlet convolutions, we see in Fig. [9]that for n = 
the NLO corrections are comparatively moderate for the imaginary part and larger for the real part. 

Going from n = to higher Gegenbauer indices n = 2 and n = 4, the NLO corrections become 
larger, as we see in Figs. [8] and [9] and already observed at small xb- Generically this is not unexpected, 
since the z dependent kernels (110p contain logarithms In z and lnz which enhance the endpoint 
regions of the z integration, and those endpoint regions are more prominent for higher Gegenbauer 
polynomials in the expansion (|15p . Note that according to phenomenological estimates or lattice 
calculations the coefficients a n of these polynomials are clearly smaller than ao, so that increasing 
corrections to 7i n for higher n do not affect the sum Y^ n a-nTin as much. We note that in the modified 
hard-scattering approach of Sterman et al. [14J, which goes beyond the collinear approximation used in 
the present work, the endpoint regions in z are suppressed by radiative corrections that are resummed 
into Sudakov form factors. As just discussed, we do not observe such a suppression in the fixed-order 
results analyzed here, where various positive and negative corrections compete with each other — only 
some of them related to the Sudakov factor. How the situation will be at higher orders is an important 
question, which goes beyond the scope of the present work. 

Let us now take a closer look at the /j,r dependence of the corrections. As we explained in Sect. [21 
the pure quark singlet kernel T& is independent of this scale at 0(a 2 ). According to ([TT]) the gluon 
kernel T g depends on fj? R only through (3q \n{^i 2 R / ^q PD ), which originates from graphs with gluon 
propagator corrections such as the one shown in Fig. [TOl The ^gpd dependence of this term is 
connected with the contribution proportional to (3q in the evolution kernel V" for the gluon GPD, 
given in (f92l) . As already pointed out in Sect. HI the term (3q ln(/j, 2 R / [1q PD ) does not contribute to the 
large- 1 ?/ 1 behavior of Imt 5jn (y) and is hence not relevant for the huge NLO corrections at small xb- 

For the kernel T a the situation is more involved. The general structure of its convolution with the 
quark singlet distribution H s can be written as 

< ^ = 00 (<« + <£> In pj + <g + <g In -f- + <5 In (59) 

V /*R/ ^GPD ^DA 

(3) (8) 

with an analogous decomposition for the convolutions Tin and Tt n . The terms proportional to (3q 
originate from graphs with gluon propagator corrections such as in Fig. [TU1 whereas the terms with 
subscripts C, G, D do not contain /3q. In Fig. [TT] we show the corresponding contributions for n = 0. 
We see that terms multiplying ln(Q 2 / ^q PD ) are rather small, whereas those going with ln(Q 2 / fi 2 DA ) 
are of course absent for n = 0. The term H.o,R is clearly smaller than TLq^ and has opposite sign. 
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Figure 8: LO and NLO terms of the convolutions in the gluon and quark singlet sector at Q = 2 GeV, 
with n R = fi GPD = hda = Q- 
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Figure 9: LO terms and the sum of LO and NLO terms of the convolutions in the quark non-singlet 
sector at Q = 2GeV, with fj,R = ^gpd = fJ-DA = Q- 




Figure 10: Example graphs giving rise to terms proportional to [3q in the hard-scattering kernels T a 
and T g . 
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Figure 11: Individual terms (|59p in the convolutions of the quark non-singlet kernel for n = 0. 
The quark distributions are evaluated at hgpd = 2GeV and the running coupling in the kernels at 
fi R = 2GeV. 

TLo^c has also the opposite sign compared to TCq^ but is similar in magnitude. We note that \TC n ./3\ 
increases with n, as can be seen from (j86|) and (|88p . 

Let us briefly comment on the BLM scale setting prescription [23], which has been discussed in the 
context of exclusive meson production in \1Q\ [24"] . This prescription aims at including the corrections 
from graphs like those of Fig.^Ojin the argument of the running coupling, and for the case at hand takes 
Hr such that the contribution from TL n ,i3 cancels against the one from Tt n> R ln(Q 2 /fj, R ) in ([59]) , As is 
evident from Fig. [XT] this requires fj, R to be substantially lower than Q 2 . For most of experimentally 
accessible kinematics, the resulting fj,R is in fact far below the region where perturbation theory 
can be applied. In such a situation, the perturbative running of a s is often modified such that the 
coupling saturates for decreasing fiR. We note that in the context of our NLO analysis, the logarithm 
00 ln(Q 2 / fi R ) in the hard-scattering kernel is intimately related with the perturbative running of 
a s (nR), so that keeping one while modifying the other is not obviously consistent. 

We also remark that if TL n ^ and TL n ,R m (Q 2 /V|j) are ma de to cancel by the BLM scale choice, 
one is left with a relatively large correction from TL n) c- For scale choices where \x R is closer to Q 2 , one 
instead has a partial cancellation between TC n; c and H. n p. A more detailed analysis for the similar 
case of the electromagnetic form factor is given in [25J, which also discusses the issue of Sudakov-type 
corrections we raised above. 

Figure [12] shows the dependence of the convolutions on fiR at fixed hgpd = Q- Within the 
range shown we generally find a moderate scale dependence, both at LO and at NLO. An exception 
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Figure 12: Dependence of the convolutions on /i^j. The factorization scale is held fixed at \xqpd = Q- 
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Figure 13: Dependence of the convolutions on hgpd- The renormalization scale is held fixed at 
m = Q- 
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Figure 14: Dependence of the convolutions on the common scale fj, = fJ,R = ^gpd = Hda- 
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Oi(fiDA)/a 2 {no) x Re (W 9 + W s(o) + H s{b) ), x B = 0.18 a 2 (n DA )/a 2 (n ) x Im (ft? + H s{a) + H s(6 >), x B = 0.18 
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Figure 15: Dependence on the factorization scale fin A of the convolutions in the gluon and quark 
singlet sector multiplied by the scale dependence (^da) / 0-2(^0) of the corresponding Gegenbauer 
coefficient. The reference scale for 02 is taken as no = 2 GeV, and the other scales are set to 
[i-GPD = VR = Q = 2 GeV. 

is the region /^#<2GeV , where the growth of the LO results simply reflects the growth of a s (fiR)- 
Note that with the parameters specified at beginning of Sect. H] we have a^ 3 ^(2GeV) = 0.30 and 
c4 3 )(lGeV) = 0.51. The NLO results further contain explicit logarithms ln(Q 2 //^), which in some 
cases can cause corrections to grow out of control, especially for the real parts of convolutions. We 
note that for HeTi^ the NLO correction is is unusually large compared with the LO term. This is 
because of a nearby zero in ig, as is seen in Fig. [9] and should not be a reason of particular concern. 

The variation of the convolutions with hgpd at fixed fiR = Q is shown in Fig. [13] We again 
find a rather moderate scale dependence, except when [Igpd becomes small. The dependence on a 
single scale fi = fj,R = n-GPD is shown in Fig. [T41 Note that in many cases the individual variation 
of fj,R decreases the amplitude in absolute size whereas the variation of ^gpd increases it, with both 
tendencies partially canceling when the scales are set equal. Again we find that the scale dependence 
becomes quite drastic below 2 GeV. 

We finally discuss the dependence on fir)A for Gegenbauer indices n > 0. According to (fT8|) 
the convolutions TCnil^DA) appear multiplied by a n (/iDA) in the process amplitude, where the scale 
dependence of both factors partially cancels. In Fig. [15] we therefore plot convolutions multiplied 
with a.2(fiDA)/ 02(^0) = [ct s (fJ'DA)/o: s (fio)]' y2 ^ following the relation (fTBT) , The corresponding plots 
for n = 4 and for convolutions in the quark non-singlet sector look very similar. We find that the 
dependence on [ir>A is slightly decreased when going to NLO. 

6 Proton helicity flip amplitudes 

We now turn to the convolutions of the hard-scattering kernels with the GPDs describing proton 
helicity flip. In this section we take t = —0.4 GeV 2 , which is the value for which will present estimates 
for observables in the next section. 

In Figs. [16] and [17] we compare the LO and NLO terms of the convolutions in the gluon and quark 
singlet sector for the two models described in Sect. 13.11 For model 1 the individual corrections for 
gluon and quark convolutions look quite similar to those we saw for 7i in the previous section. The 
sum of gluon and quark singlet contributions at LO is however very small in this model because 
of cancellations, so that the NLO term dominates in a wide kinematical region. In model 2 the 
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Figure 16: LO and NLO terms of the convolutions in the gluon and quark singlet sector for model 1 
at Q = 2GeV and t = — 0.4GeV 2 . The scales are set to fj,R = /j-cpd = Q- 
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Figure 18: LO and NLO terms of the convolutions in the quark non-singlet sector for model 1 at 
Q = 2 GeV and t = -0.4 GeV 2 , with fx R = fj, GPD = Q. 
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Figure 19: As Fig. [T8l but for model 2. 
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gluon contributions are nearly absent, so that the quark singlet contribution dominates in this sector. 
We note that, contrary to the individual terms, the sum of gluon and quark singlet contributions 
comes out to be rather similar in the two models and is small compared with the flavor non-singlet 
contributions shown in Figs. [TH] and [191 According to our discussion at the end of Sect. 13.11 this has 
its origin in the sum rule (|43ft for the second moment of E at t = 0, so that we expect a small net 
contribution from gluons and the quark singlet in large class of models for E. 

As shown in Figs. [18] and [I9j the NLO corrections to the quark non-singlet convolutions are 
relatively moderate but not small, similarly to the case of Ti. The size of the convolutions is quite 
different in the two models, indicating the important role played at intermediate xb by sea quarks in 
model 1. Let us recall that with the double distribution ansatz ([33]) the GPDs at x ~ £ are sensitive 
to forward parton distributions with momentum fractions well below £, as discussed in Sect. 4.3.3 

of HZ]. 



7 Cross sections and asymmetries 

Having studied in detail the building blocks of the scattering amplitude for vector meson production, 
we now combine them to observables. We recall that to leading order in 1/Q there are just two of 
these: the unpolarized 7*p cross section and the asymmetry for a transversely polarized target, both 
referring to longitudinal polarization of virtual photon and produced meson. The ep cross section in 
the leading 1/Q approximation can be written as 



da(ep — > epV) a 1 
dtdQ 2 dyd<fid<t)s ~ 47r 3 



Q' 2 



xb 1-y dai 
V dt 



1 + St sin(0 - <fis) A UT 



(60) 



where y is the usual inelasticity variable for deep inelastic scattering and St denotes the transverse 
component of the target polarization, (j) is the azimuthal angle between lepton plane and hadron plane, 
and 4>s is the azimuthal angle between lepton plane and target spin vector, both defined according 
to the Trento convention [40] . The ^*p cross section da^/dt and the polarization asymmetry Aut 
depend on xb, Q 2 and t. To leading order in 1/Q they are given by 



d(JL 

~~dT 



7T 



a 



XB 



1- X B 



fv 



(1 



f)\n v \ 2 



(t/(4m 2 ) +£ 2 ) \£ v \ 2 - 2£ 2 Re(^7A 



(61) 



and 
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VW 2 lm{£^n x 



where to 



(1 - £ 2 ) \H V | 2 - (t/(4m|) + e) \£ v \ 2 - 2£ 2 Re(£ v H v ) ' 
-4m,p£ 2 /(l — £ 2 ). Here we have combined the convolutions (|19p into 

OO i- 



(62) 



n=0 



(63) 



with analogous combinations for Tiy and £y . In the remainder of this section we take the asymptotic 
form of the meson distribution amplitude, i.e. we set a n = for n > 2. As long as £y is not much 
larger than Tiy, the cross section (|6ip is dominated by the term with |7Yv | 2 in a wide range of 
kinematics, where the prefactors £ 2 and i/(4m 2 ) of the other terms are small. The asymmetry (I62D 
is then approximately given by 
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sin 5y 



(64) 
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Figure 20: Real and imaginary parts of the convolutions (|63p for p production at Q = 2 GeV and 
t = —0.4 GeV 2 , with model 1 taken for the proton helicity-flip distributions. The scales are set to 
= hgpd = Q- 



K, at Q 2 = 4GeV 2 t = -0.4 GeV 2 

Re LO 

Re LO+NLO 

Im LO 

Im LO+NLO 



0.1 



0.2 



0.3 



0.4 



1 

-1 
-2 
-3 
-4 
-5 
-6 



£ w at Q 2 = 4 GeV 2 , t = -0.4 GeV 2 




0.1 



0.2 



xb 



0.3 



Re LO 

Re LO+NLO 
Im LO 

Im LO+NLO 

' oa' 



Figure 21: As Fig. [5D] but for lo production. 



where 5y = aigiTLy / 8y) is the relative phase between TLy and Ey. 

Figure [2D] shows the real and imaginary parts of the convolutions appearing in (|61|) and (|62p . 7i p 
is dominated by the gluon and quark singlet part, and in line with our discussion in Sect. [5] we find 
rather moderate corrections for the imaginary part but a very unstable real part in a wide range of xb ■ 
As for £ p , its real part is very small and subject to large relative corrections, whereas its imaginary 
part is much larger and receives corrections of order 100% . As we see in Fig. [T8l the individual flavor 
non-singlet combinations 8^ and 8^ are less affected by corrections, but they have opposite sign 
and partially cancel in the sum relevant for p production. The rather small but unstable contribution 
from the gluon and quark singlet terms is hence important in this channel and largely responsible for 
the NLO corrections seen in Fig. [20l As a further consequence of the cancellations just mentioned, 
the size of 8 P is tiny compared with TC P . The quark flavor combination relevant for p production is 
2u + d, where in our model the flavor combinations add for H but largely cancel for E. 

For uj production we see in Fig. [21] that at smaller xb the convolution Ti^ is about 1/3 of TC P , 
which follows from the dominance of the gluon contribution as seen in ([6]) and ([7]), whereas at larger 
xb differences between the two channels appear. £ w is much bigger than 8 P since in the combination 
2u — d the contributions of u and d quarks add for E, and the size of its radiative corrections reflects 
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the one of 8^ in Fig. [THl We note that the dominance of the imaginary over the real part in and 
£ p is less pronounced in model 2, as can be anticipated by comparing Figs. flBI and [I~9l 

The cross section dai/dt for p production is dominated by (ImWp) 2 , except for contributions 
from (ReTtp) 2 at small xb for NLO and at large xb for LO. Given the size of corrections to lmTi p 
in Fig. [20] we thus have quite substantial NLO effects in the cross section at Q 2 = 4GeV 2 , as shown 
in Fig. [22] For Q 2 = 9GeV 2 and xb > 0.1 the relative corrections decrease. The plot has been 
calculated with model 1 for E, but since its contribution to dai/dt is negligible the corresponding 
curves for model 2 look very similar. To obtain an estimate of scale uncertainties, we show bands 
corresponding to fj, = [ir = p-gpd between 2GeV and 2Q. Given our discussion in the previous 
section, we do not consider it meaningful to go to scales below 2 GeV, so that the bands in the figure 
are strongly asymmetric. For Q 2 = 4 GeV 2 they go only in one direction, and the band of the LO 
result does not provide an estimate for the size of the NLO corrections, which turn out to go in the 
other direction. 

We have a very peculiar situation for the polarization asymmetry Ajjt in P production, which 
as shown in Figs. [23] and [Ml is very small in both models 1 and 2 due to the cancellations in £ p 
discussed above. Ajjt changes quite dramatically from LO to NLO in a wide range of kinematics, 
clearly because of the NLO corrections in the numerator. A closer look at Fig. [20] reveals that the 
large perturbative corrections in Im(£* Tip) are mainly due to the large corrections to both KeTL p 
and Re£p. These hardly affect the unpolarized cross section, which is strongly dominated by ImW p . 
At higher Q 2 the instability of Ajjt is less pronounced, and in model 2 we even have quite small 
corrections. We note that the bands from the scale variation at LO order are extremely narrow in 
Figs. [23] and [23] This is because the scale variation of a s {pLfi) cancels in the ratio Ajjt at LO and 
because in the kinematics we are looking at, the pgpd dependence of 7i p and £ p is rather weak. In 
this situation, the scale uncertainty of the LO result does obviously not provide a good estimate for 
the size of higher-order corrections. Let us finally remark that at t = —0.4 GeV 2 the asymmetry Ajjt 
must go to zero as tends to 0.484 because of the prefactor \fto — t in ([62]) . 

The cross section for uj production is shown in Fig. [25] and shows a similar pattern of NLO 
corrections to the one in p production, reflecting the similar pattern of corrections we have seen for 
ImWp and Im^. As a result the ratio of cross sections da^jdt in the two channels is quite stable 
under radiative corrections, as seen in Fig. [26] The target polarization asymmetry, shown in Fig. [27] 
for model 1, changes however drastically between LO and NLO at small to intermediate xb- This is 
because ImS^ then dominates over R,e £^j , so that its product with the unstable convolution R,e 
controls the numerator of the asymmetry. The absolute size of Ajjt can be large in this channel 
since \8 U \ ~ [H w \ in our model. According to Fig. [2H the relative phase 5 U is close to zero at LO for 
xb<0.3, so that the factor sin^ in (|64p makes Ajjt small and prone to large radiative corrections. 

Let us finally take a look at <fi production. At LO this channel is strongly dominated by gluon 
exchange, since in our models strange quark distributions are small for H and even more so for E. 
At NLO we have further contributions from the pure sing let terms H s( ~ b) and £ s(fe) , which are not 
negligible. We see in Fig. [28] that the NLO corrections to the cross section are large at small xb and 
slowly decrease with xb- Except for the region of small xb, this pattern is quite different from the 
one in p production, so that the cross section ratio for the two channels receives important corrections 
at larger xb as we see in Fig. [29] The asymmetry Ajjt is essentially zero at LO, because in our model 
the relative phase 6a between Ti^ and £§ is very close to zero. This changes at NLO, where in model 1 
we obtain a small to moderate Ajjt, as shown in Fig. [30] 
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Figure 22: Longitudinal cross section for 7*p — » in model 1. Bands correspond to the range 
2 GeV < \i < 4 GeV in the left and to 2 GeV < \i < 6 GeV in the right plot, and solid lines to fi = Q 
in both cases. 
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Figure 23: The transverse target spin asymmetry defined in (|60p . calculated for model 1. The meaning 
of the bands and solid lines is as in Fig. [22j 
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Figure 24: As Fig. but for model 2. 
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Figure 25: As Fig. [22] but for j*p — * up. 
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Figure 26: Ratio of cross sections da^jdt for u) and p production in model 1. The meaning of the 
bands and solid lines is as in Fig. | 
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Figure 27: As Fig. [23] but for -y*P ~^ ^p. 



37 



da L /dt [nb/GeV 2 ] at Q 2 = 4GeV 2 , t = -0.4GeV ; 
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Figure 28: As Fig. [22] but for Yp ~^ 4>P- 
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Figure 29: Ratio of cross sections daL/dt for <p and p production in model 1. The meaning of the 
bands and solid lines is as in Fig. | 
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Figure 30: As Fig. [23] but for Yp ~^ 4>P- 
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8 Pseudoscalar meson production 



Having studied in detail the production of vector mesons, let us finally take a look at pseudoscalar 
production. We will only consider ^y*p — > 7r + n, which was already studied at NLO in [10J. Gluon 
distributions do not contribute in this channel. 

In the collinear approximation the amplitude for this process can be written as 
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(65) 



n=0 



with e u = 2/3, = —1/3 and f n = 131 MeV. <f>ir(z) is the twist-two distribution amplitude of the 
pion and has a Gegenbauer decomposition as in (fT5|) . The convolutions T£ are defined as 



dx 



e u T a , n (x, - e d T a (-x, £) F u (x, £, t) - F d (x, £, t) 



(66) 



and the kernels T a (z,x, £) and T ajra (x,£) are the same as in Sect. [2j The matrix elements F q are the 
counterparts of F q for polarized quarks and given by 
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2m r , 



(67) 



in terms of the generalized parton distributions H and F, where as in the unpolarized case we use 
the conventions of [T7]. Since the hard-scattering kernel in (|66p is neither even nor odd in x, the 
convolution involves both the charge-conjugation even and odd combinations 



F*« (x, C t) = F*(x, e, t) + F*(-x, e, t) , F'H (x, t) = F*(x, £, t) - F«(-x, £, t) 
We model the distributions H in close analogy to the unpolarized case and set 

"1-1/31 
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(69) 



with /i( 2 )(/3, a) as in (|24p and 

ff«M(x,0,i) = A^(x)exp[i/ ?I ,(x)] +2Ag(x)exp[t/ 9 - 
ff g (-)(x,0,t) = Aq v (x)exp[tf qv (x)] 



(70) 



for x > 0. The values for x < are determined by the symmetry properties following from (|68|) . For 
the polarized valence and antiquark densities Aq v and Aq we use the NLO parameterization from 
[4T] at fj, = 2GeV, and for the t dependence we take the same functions fq v (x) as in ([26]) . ([2?]) and 
furthermore set fq{x) = fq v {x). As was shown in [28J, this gives a good description of the isovector 
axial form factor via the sum rule 



F A (t) = f 1 dx (a, 0,t) - H d (+\x,0,t)} . 

Jo 



(71) 
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0.2 0.3 



Figure 31: The convolution TLq denned as in (1661) . evaluated at Q = 2GeV and t = -0.4 GeV 2 . The 
scales are set to /j,r = [igpd = Q- 



For the nucleon helicity-flip distribution E we take a pion exchange ansatz 



E u (x,a,t) = -E d (x,a,t) 



»(g-M) 

2? 



x + £ \ ^ P 9A A 2 
2£ ) ml 



mZ 



t A 2 -t ' 



(72) 



with the nucleon axial charge gA ~ 1-26 and a cutoff parameter A = 0.8 GeV [42j to suppress large 
off-shellness of the exchanged pion in the t channel. 



8.1 Results 

The convolution at LO and NLO is shown in Fig. [31] for n = 0. We find moderate corrections for 
the imaginary part and larger ones for the real part. For we can easily give the analytic form of 
the NLO result. The scale dependent terms admit a closed expression, 



Q n £ n OC ^ a m a n \ 1 + ^ 

n rn n >■ 



R , 14 7m + 7n 



In 



/4 



7m In 



Q 2 
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7nln 



+ ... 



(73) 



where the . . . denote contributions which depend neither on Q 2 and the scales nor on (3q. Including 
these terms we can write 



^a n £^ oc (1 + a 2 + a 4 ) 2 + 



7T 
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(74) 



where we have set nj = 3 in [3q and where we approximated numerically the coefficients written 
with a decimal point. Here the . . . denote terms with higher Gegenbauer coefficients. Note that 
these coefficients appear twice, once for the produced pion and once for the pion exchange ansatz 
of the distribution E. Up to a global factor, the expression (|74|) also gives the NLO result for the 
electromagnetic pion form factor F n (Q 2 ) at large spacelike momentum transfer Q 2 , and it agrees with 
the result in the detailed study [9] . Let us first discuss the case m = n = relevant for the asymptotic 
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form of the pion distribution amplitude, where the convolution has no dependence on hgpd and /j,da- 
We then have the rather large coefficient 79/12 ~ 6.6 in square brackets, so that with the scale choice 

= Q there are quite large NLO corrections. The corrections are zero for ii 2 R = e _79 / 27 Q 2 ~ 0.05 Q 2 , 
which is outside the perturbative region for most cases relevant in practice. The BLM scale for this 
case is yet smaller: with ([73]) we reproduce the well-known result = e _14//3 Q 2 « 0.01 Q 2 [10\ . The 
coefficient of a s j-K in ([74]) is then —47/12 —3.9 and thus again rather large, but of course the scale 
fj,R is outside the perturbative region for all experimentally relevant kinematics. We finally see in (I74h 
that for higher Gegenbauer moments the correction terms are larger than for m = n = 0. The reason 
for this is the same which we discussed in Sect. [5] for the convolutions 7i. In (j73|) we also see that the 
BLM scale becomes smaller for higher m and n. 

The observables for exclusive pion production at leading order in 1/Q are the same as for vector 
meson production, and the ep cross section is given as in (]60p . The cross section for a longitudinal 
photon and the transverse target asymmetry are now respectively given by 



dt 



7T 
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(2- XB y 

1-XB 
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with 



(77) 



n=0 



n=0 



For numerical estimates we take the asymptotic pion distribution amplitude in the following, setting 
a n = for n > 2. We note that the recent lattice study [22] obtained a rather moderate value 
a^o) = 0.201(114) at fiQ = 2GeV. 

In Fig. [32] we show the separate contributions from the terms with \Tt w \ 2 and with \£ n \ 2 in (|75p . 
as well as the full result. We see that at the value of t chosen here, the contribution from \H-k\ 2 is 
more important, mainly because of the suppression factor (A 2 — m 2 )/(A 2 — t) in our model (|72p for 
E. The square of this factor is 0.36 at t = — 0.4GeV 2 . 

We compare the LO and NLO results for the cross section in Fig. [33] and find that the NLO 
corrections are quite large, even at Q 2 = 9GeV 2 . In contrast, the corrections for the beam spin 
asymmetry are very small as seen in Fig. I34"[ in line with the findings reported in [10]. Note that with 
our model £ n is purely real, so that at intermediate xb the large relative NLO corrections in He7i. n 
do not affect the numerator of Aut in fl75|> . Approximating the asymmetry as 



A 



UT 



YHt. 



sin^T,- 



(78) 



with 5 n = arg(TC n /£ n ), we can understand why only small corrections are seen in this case: the 
relative phase 5^ is well different from zero, and the NLO corrections increase both \TCtt\ and \£n\- 
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Figure 32: The longitudinal cross section for ^y*p — ► 7r + n, evaluated at NLO. Shown are the separate 
contributions from the terms with l?^ 2 and with \£ n \ 2 in (fT5j) . as well as the complete expression. 




Figure 33: Longitudinal cross section for "f*p — ► 7r + n. Bands correspond to the range 2GeV < \i < 
4 GeV in the left and to 2 GeV < fi < 6 GeV in the right plot, and solid lines to \x = Q in both cases. 



Ajjt at Q 2 = 4 GeV 2 , i = -0.4 GeV 2 A UT at Q 2 = 9 GeV 2 , t = -0.4 GeV 2 




Figure 34: The transverse target spin asymmetry for tt + production, as defined in (|60p . The meaning 
of the bands and solid lines is as in Fig. [33j 
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9 Summary 



In this work we have analyzed the NLO corrections for exclusive meson production at large Q 2 in the 
collinear factorization approach. Using the Gegenbauer expansion of meson distribution amplitudes, 
we have rewritten the hard-scattering kernels of [11] into functions depending on only one variable, 
and we have separated the explicit logarithms in the factorization scale for the meson distribution 
amplitude and the generalized parton distributions. 

For vector meson production at small ig we find huge NLO corrections even for Q 2 well above 
lOGeV 2 , in agreement with the results obtained in [11] . The corrections have opposite sign compared 
to the Born term and can be traced back to BFKL type logarithms in the hard-scattering kernels, 
which appear with rather large numerical prefactors in this process. We conclude at this stage 
that a quantitative control of radiative corrections at small xb will require resummation of these 
logarithms. First steps in this direction have been reported in [43J. If successful, such a resummation 
in combination with a dispersion relation [39J may also be useful for stabilizing the real part of the 
amplitude, where we find very large NLO corrections even at xb ~ 0.1. 

At intermediate to large xb, typical of fixed-target experiments, we have investigated the pro- 
duction of p , oj, 4> and of tt + . We find NLO corrections to the longitudinal cross sections of up to 
100%, which somewhat decrease in size when going from Q 2 = 4GeV 2 to 9GeV 2 . Note that the me- 
son production cross section depends quadratically on generalized parton distributions — the increased 
sensitivity to these basic quantities comes with an increased sensitivity to higher-order corrections. 
We generally find that uncertainties on the cross section due to the choice of renormalization and 
factorization scales are not too large at LO and do not significantly decrease when going to NLO. For 
scales below 4GeV 2 , however, NLO corrections often grow out of control. The cross section ratio for 
to to p production turns out to be very stable under corrections, but less so the one for <p to p. For 
the transverse target polarization asymmetry Ajjt in 7r+ production we find quite small NLO effects, 
confirming the results in [10]. For vector meson production this is however not the case. With the 
models we have used for the nucleon helicity-flip distributions E, the numerator of the asymmetry in 
this channel is dominated by the product (Im£y)(ReHy) in a wide range of kinematics and therefore 
suffers from the perturbative instability we find for HeTiy at small to intermediate xb, even if the 
corrections to Im£y are not too large. It is often assumed that corrections tend to cancel in asym- 
metries. The examples we have studied show that this may hold in specific cases but not in others, 
and that special care is needed for observables like Ajjt that depend on the relative phase between 
amplitudes. 

We should recall that in the kinematics we studied, one must expect that our leading-twist results 
receive power corrections that cannot be neglected when comparing with data. They will certainly 
affect the cross sections and will not always cancel in cross section ratios. An example is the transverse 
target polarization asymmetry in ir + production. The phenomenological estimates in [12] found that 
the convolution 7i n is decreased by effects of transverse parton momentum in the hard scattering, 
whereas £ n is increased by the soft overlap mechanism that has been extensively studied in the context 
of the pion form factor. Together, these corrections may significantly increase leading- twist estimates 
for Aut- 

From our numerical studies we must conclude that a precise quantitative interpretation of exclu- 
sive meson production requires large Q 2 , say above lOGeV 2 . In addition it would be highly valuable 
to have a consistent scheme for combining radiative with power corrections, at least in parts. Nev- 
ertheless, we find that valuable information on generalized parton distributions can be obtained also 
from data at lower Q 2 . In particular, a large measured asymmetry Ajjt in vector meson production 
would give valuable constraints on the size of the proton helicity-flip distribution E 9 for gluons, which 
are most difficult to obtain in deeply virtual Compton scattering or from lattice QCD calculations. 



43 



Acknowledgments 



We gratefully acknowledge discussions with L. Favart, H. Fischer, P. Kroll, A. Rostomyan and 
A. Schafer. Special thanks are due to D. Yu. Ivanov for numerous conversations and advice. This 
work is supported by the Helmholtz Association, contract number VH-NG-004. 



A Polylogarithms 

We collect here some properties of the polylogarithms that appear in the hard-scattering kernels for 
meson production. Their definitions are 



Li 2 z 
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Li 3 z = / — Li 2 {zt) , 
Jo * 



from which one readily obtains for the imaginary parts 
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The limiting behavior for y — > —00 can be obtained from the expansions 
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which are valid for y < — 1, and from 
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which holds for y < 0. A useful relation finally is 
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A wealth of further information can be found in 



(83) 



B Hard-scattering kernels for higher Gegenbauer moments 



In this appendix we give the analogs of the hard-scattering kernels in (|2ip for Gegenbauer index n = 2 
and n = 4. For the gluon kernel we find 
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and for the pure singlet kernel 
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The quark non-singlet kernel reads 
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of the anomalous dimensions, we can give a closed form for the scale dependent terms for all even n, 
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where the terms denoted by . . . are independent of Q 2 and the scales and do not involve 0q. From 
the scale dependence f|16j) of the Gegenbauer coefficients of the meson distribution amplitude we can 
readily reconstruct their evolution equation 
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With ([20]) and ([88]) we see that the \i da dependence of the process amplitude (fTBj) cancels up to terms 
of order a 3 , as it must be. 

C Evolution kernels 



For definiteness we give here the LO evolution kernels for GPDs, which we have used to check the 
scale invariance of the NLO amplitude for meson production as explained in Sect. [21 The non-singlet 
evolution equation reads 



^ pNS{x > *> t} = £ I yNS (f * f ) pNS{y > *> t} • 
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where F NS can be a flavor non-singlet combination such as F u ^ — F d ( + \ or the charge-conjugation 
odd combination F q( -\x, f, t) = F q (x,£,t) + F q (-x,£,t) for a single quark flavor. In the gluon and 
quark singlet sector we have a matrix equation 
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with Tp = 1/2 and the remaining constants as given in (|13|) . The plus-prescription appearing in V qq 
and V ss is defined by 
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and the function p(x, y) specifies the support as 
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The evolution equations for polarized GPDs read as in (|90p and (|9ip . with the unpolarized matrix 
elements F and kernels V replaced by their polarized counterparts F and V. With F q ( + *> and F q<y ^ 
defined in flSSJ above, F^ 5 can be either a flavor non-singlet combination like F u( +) — or a 

charge-conjugation odd combination F q (~\ whereas the flavor singlet combination is given by 



F s _ p u{+) + p d{+) + ps{+) _ 

To 0(a s ) the polarized evolution kernels are 

V NS (x,y)=V qq (x,y) = V qq (x,y) 
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The kernels given here agree with those in [45] if one takes into account that any contribution to 
V 9q (x,y) which is even in y at fixed x will drop out in the convolution (|9ip . Taking the limit £ — > 

as 

ijif Fw (t'0 ^(*.m_pw (98) 

one obtains the usual DGLAP evolution kernels from (|92j) . and in analogy one recovers the polarized 
DGLAP kernels from ([97]) . The factors j in front of V s9 and V" reflect the different forward limits 
of the quark and gluon GPDs. 
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